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SECTION  I 


INTRODUCTION 


This  document  represents  the  final  report  of  the  Satellite  Image  Quality 
Simulation  Study,  performed  for  Rome  Air  Development  Center  under  Contract 
No.  F30602-74-C-0209. 

1.1  PURPOSE 

The  purpose  of  the  study  was  to  develop  a quantitative  model  for  the  re- 
lationship between  the  quality  of  images  of  satellites,  obtained  through  a 
turbulent  atmosphere  from  large  ground-bas«_d  telescopes,  as  perceived  by  human 
observers,  and  quantitative  measures  of  thi  optical  transfer  function  (OTF) 
and  noise  characteristics  of  the  image  formation  process.  To  determine  this 
relationship,  a series  of  images  of  satellites  was  generated  in  which  the 
effects  of  various  OTF  and  noise  degradations  were  simulated.  These  simula- 
tions were  produced  by  first  digitizing  a high-resolution  undergraded  photo- 
graph of  a satellite,  simulating  noise  and  OTF  degradations  in  a digital 
computer,  and  then  reconstructing  the  degraded  images  on  a high-quality  image 
playback  device. 

The  images  were  then  viewed  by  a group  of  human  observers,  who  were 
asked  to  compare  the  degraded  satellite  images,  one  at  a time,  to  a set  of 
reference  images  of  the  same  satellite  having  known  resolution  characteris- 
tics; after  each  simulated  image  was  viewed  by  the  group  of  observers,  the 
relative  quality  of  that  image  (relative  to  the  reference  set)  was  statis- 
tically determined  by  computing  the  average  quality  estimate  from  all  ob- 
servations. 

The  reference  set  consisted  of  a series  of  images  in  which  the  diffrac- 
tion effects  of  various  aperture  diameters  ranging  from  60  inches  to  4 inches 
were  simulated.  The  observed  quality  of  the  simulations,  as  determined  by 
the  observers,  is  expressed  in  terms  of  the  equivalent  aperture  diameter  of 
the  member  of  the  reference  set  to  which  the  observed  image  was  equated. 
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The  results  of  these  viewing  and  interpretation  (V  and  I)  experiments 
were  then  compared  to  the  predictions  of  several  image  quality  models,  and 
the  model  that  produced  the  best  fit  to  the  observed  data  was  selected. 

1.2  SUMMARY  OF  EXPERIMENTS 

The  main  goal  of  the  program  was  to  generate  images  for  the  purpose  of 
fitting  the  results  to  a model,  and  not  to  simulate  the  output  of  any  specific 
hardware  system;  therefore,  the  images  were  generated  with  a wide  range  of 
noise  and  OTF  degradations,  solely  to  provide  as  large  a base  of  images  as 
possible.  However,  the  ultimate  goal  of  the  Teal  Blue  Program  was  kept  in 
mind  - i.e.,  the  fabrication  of  a Compensated  Imaging  System  in  which  both 
pre-compensation  of  the  wavefront  is  performed  before  the  images  are  recorded, 
and  post-compensation  is  performed  on  the  images  after  they  are  recorded. 

With  this  operational  concept  in  mind,  the  noise  mechanisms  and  the  OTF 
degradations  selected  for  the  simulations  were  designei  to  be  the  generic 
types  of  degradations  that  a compensated  imaging  system  might  suffer.  For  the 
noise  mechanism,  both  multiplicative  noise  and  additiv*  noise  were  simulated; 
for  the  OTF  degradations,  two  types  of  transfer  functions  were  generated,  one 
for  short  time  exposure  conditions,  the  other  for  long  tine  exposure  condi- 
tions. The  relative  strengths  of  the  noise  and  the  OTF  degradations  were 
varied,  and  the  operation  of  two  types  of  post-processing  restoration  opera- 
tions were  also  included  in  the  simulations.  The  first  type  of  restoration 
filter  that  was  simulated  assumed  complete  knowledge  of  the  degrading  atmos- 
pheric OTF;  for  the  case  of  short  time  exposure  conditions,  this  represents 
a situation  that  is  probably  unachievable  in  the  real  world,  but  it  is  useful 
in  placing  an  upper  bound  on  the  ability  of  post-compensation  to  perform  its 
function,  in  addition  to  providing  images  against  which  various  image  quality 
models  can  be  matched. 

The  second  type  of  post-compensation  processing  that  was  simulated  was 
the  use  of  a post-compensation  operator  that  represents  only  partial  know- 
ledge of  the  degradation  transfer  function;  this  filter  was  derived  on  only 
a statistical  knowledge  of  the  degrading  OTF,  and  assumed  that  the  degrading 
OTF  was  a smooth  transfer  function  with  no  phase  component.  Thus,  the  post- 
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compensated  images  produced  by  this  filter  contained  the  effecte  of  uncor- 
rected modulus  and  phase  errors  in  the  OTF. 

Figure  1 illustrates  the  various  combinations  of  the  degradations  men- 
tioned above  that  were  used  to  generate  the  simulations;  Table  1 lists  the 
number  of  images  generated  during  the  study.  A total  of  121  images  were 
generated  and  used  in  the  study. 

TABLE  1.  SUMMARY  OF  EXPERIMENTS 

Set  Experiment 

A Various  Diffraction-Limited  Apertures 

B Diffraction-Limited  Apertures  and  Poisson  Noise 

H Diffraction-Limited  Apertures  and  Gaussian  Noise 

C Short  Time  Exposure  OTF  and  Poisson  Noise 

I Short  Time  Exposure  OTF  and  Gaussian  Noise 

D Short  Time  Exposure  OTF  and  Poisson  Noise  and  Post- 

Processing  (Exact) 

K Short  Time  Exposure  OTF  and  Gaussian  Noise  and  Post 

Processing  (Exact) 

J Short  Time  Exposure  OTF  and  Poisson  Noise  and  Post- 

processing (Partial) 

F Long  Time  Exposure  OTF  and  Poisson  Noise 

G Long  Tine  Exposure  OTF  and  Poisson  Noise  and  Post- 

Processing  (Exact) 


No . of 
Images 

15 

12 

12 

24 

9 

9 


J 
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Figure  1.  Sunniary  of  Experiments 
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SECTION  II 
IMAGE  SIMULATION 


In  this  section,  the  operations  performed  in  the  generation  of  the  images 
used  to  perform  the  viewing  and  interpretation  experiments  will  be  discussed. 

2.1  IMAGE  GENERATION 

The  source  of  all  of  the  images  generated  during  the  study  was  a photo- 
graph of  a model  of  the  P72-2  satellite  that  was  taken  in  a solar  simulator; 
the  model  was  produced  by  Rockwell  International  Corporation,  under  Contract 
No.  F30602-74-C-011. . 

The  original  negative  of  the  satellite  model  was  used  to  produce  a large- 
format  positive  transparency  suitable  for  scanning  on  Perkin-Elmer ' s Line  Scan 
Image  Generator  (LSIG) , which  was  used  to  digitize  the  image. 


2.1.1  Scanning  Parameters 

The  positive  transparency  was  scanned  at  a spatial  rate  that  corresponded 
to  the  sampling  of  a 6.5-arc-second  satellite  at  the  Nyquist  rate  for  a 60- 
inch  (1.524  meter)  telescope  operating  at  a wavelength  of  0.5  micrometer.  The 
Nyquist  angular  sampling  frequency  is  given  by 


n 

s 


2D 

X 


(samples/radian) 


(1) 


which  corresponds  to  an  angular  subtense  for  each  image  picture  element 
(pixel)  of 


0 

s 


_X_ 

2D 


(radians/ sample) 


(2) 


For  the  parameters  listed  above,  each  pixel  subtends  1.64  x 10  radians,  or 
0.034  arc-second.  For  a 6.5-arc-second  object,  this  corresponds  to  a total 
number  of  191  pixels  across  the  large  dimension  of  the  object.  As  a result, 
all  of  the  simulations  generated  for  the  study  consisted  of  256 x 256  images. 
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2.1.2  Intensity  Linearization  of  the  Digitized  Image 

Because  of  the  various  nonlinear  properties  of  the  two  film  9teps  (the 
exposure  and  development  of  the  original  negative  and  the  positive  transpar- 
ency) and  the  digitization  operation,  it  was  necessary  to  linearize  the  digi- 
tized information  so  that  the  digital  values  of  the  sampled  image  were  related 
in  a linear  fashion  to  the  original  reflectivity  of  the  object  itself.  For 


tunately,  the  original  negative  contained  an  image  of  a reflective  target 
consisting  of  16  chips  of  known  reflectivities.  By  carrying  the  image  of  the 
test  chips  through  all  of  the  steps  (film  and  digitization),  it  was  a simple 
matter  to  determine  the  net  nonlinearity  of  the  combination  of  the  film  and 
digitization  operations,  and  correct  for  these  nonlinearities  by  digitally 
implementing  a look-up  table  in  the  computer  that  effectively  removed  the 
nonlinearities.  The  same  operation  was  also  used  to  optimally  scale  the 
range  of  digital  values  so  that  the  output  images  would  use  the  full  availa- 
ble dynamic  range  of  the  photographic  material  used  to  produce  the  images 
that  were  subsequently  used  for  the  viewing  and  interpretation  experiments. 


To  further  assure  consistent  quality  for  the  images  to  be  used  in  the 
experiments,  a digital  step  tablet  was  played  back  with  each  of  the  images 
generated.  This  digital  step  tablet  was  used  by  the  photo  technician  as  a 
standard  to  determine  the  optimum  exposure  for  each  of  the  negatives  produced 
on  the  LSIG  in  the  image  playback  mode.  Thus,  the  problem  of  inconsistent 
photographic  processing  introducing  errors  into  the  experiments  was  reduced 

to  a minimum. 

2.2  OPTICAL  TRANSFER  FUNCTION  SIMULATION 

For  the  study,  images  formed  through  three  types  of  optical  transfer 
functions  were  generated;  diffraction-limited  apertures  of  various  diameters, 
short  time  exposure  transfer  functions,  and  long  time  exposure  transfer  func- 
tions. The  purpose  of  simulating  the  latter  two  types  of  transfer  functions 
stems  from  the  desire  to  simulate  different  types  of  transfer  functions  for 
the  purpose  of  testing  image  quality  models  over  widely  differing  conditions. 
These  two  types  most  likely  represent  extremes  in  the  kinds  of  transfer  func- 
tions to  be  encountered  by  a Compensated  Imaging  System. 
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Short  time  exposure  images,  arc  characterized  by  the  condition  that  the 
exposure  time  for  the  image  is  short  with  respect  to  the  velocity  ot  propaga 
tion  of  the  turbulence  across  the  aperture  of  the  system.  For  this  case,  the 
turbulence  is  essentially  "frozen"  during  the  exposure  time,  and  the  result- 
ing complex  optical  transfer  function  will  in  general  have  a significant  phase 
component,  and  the  modulus  will  be  a non-monotonic  function  of  spatial  fre- 
quency and  will  have  significant  irregularities. 

If,  on  the  other  hand,  the  exposure  time  is  long  enough  to  permit  a sig- 
nificant amount  of  averaging  of  the  atmospheric  turbulence,  then  the  result- 
ant optical  transfer  function  will  have  a small  phase  component,  and  the 
modulus  of  the  transfer  function  will  be  a smooth  monotonically  decreasing 
function,  with  a lower  value  for  the  modulus  at  all  frequencies  because  of 
tne  averaging  process. 

A comparison  of  these  two  cases  for  the  condition  in  which  there  is  a 
significant  amount  of  wavefront  error  reveals  that  the  short  time  exposure 
image  is  harder  to  process  from  the  standpoint  of  restoration,  since  know- 
ledge of  the  shape  of  the  complex  OTF  at  the  instant  that  the  image  was  ob- 
tained is  needed.  On  the  other  hand,  the  long  time  exposure  image  presents 
a simpler  restoration  problem,  as  it  will  require  only  a relatively  smooth, 
real  (non-complex)  frequency  domain  filter  whose  shape  can  be  predicted  sta- 
tistically. Figure  2 contains  examples  of  these  two  types  of  transfer  func- 
tions, for  two  levels  of  mean  square  wavefront  error.  Both  of  the  top  trans- 
fer functions  have  Strehl  ratios  of  0.03,  while  both  bottom  transfer  functions 
have  Strehl  ratios  of  0.3.  These  transfer  functions,  whose  derivations  are 
discussed  in  the  following  paragraphs,  illustrate  the  differences  between  the 
two  types  of  transfer  functions.  It  should  also  be  noted  that  the  distinction 
becomes  less  distinct  as  the  overall  system  quality  improves. 

2.2.1  Reference  Set  Generation 

The  reference  set  was  generated  using  a fast  Fourier  transform  (FFT)  al- 
gorithm, for  the  series  of  15  aperture  diameters  listed  in  Table  2.  Three  of 
the  images  and  optical  transfer  functions  are  shown  in  Figure  3 and  nine  of 
the  corresponding  images  are  shown  in  Figure  B-l  in  Appendix  B.  For  the  im- 
ages in  Figure  B-l  and  all  subsequent  figures  showing  images,  the  three-digit 
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Figure  2.  Long  Time  Exposure  and  Short  Time  Exposure  OTF's 


60  Inches 


Reference  Set  Image  and  OTF  Examples 
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number  below  the  right-hand  comer  of  each  image  is  a frame  number  for  the 
image,  and  is  used  to  reference  that  image  in  all  figures  and  tables  contained 
in  this  report.  In  addition,  all  OTF  plots  shown  in  this  report  have  the  same 
scale,  and  have  as  the  maximum  frequency  displayed  the  cutoff  frequency  of  a 
60-inch  aperture. 

TABLE  2.  REFERENCE  SET  APERTURE  DIAMETERS 


Aperture 

Aperture 

Diameter 

Diameter 

,e  Number 

(Inches) 

(Meters) 

015 

60 

1.52 

018 

54 

1.37 

019 

48 

1.22 

020 

42 

1.07 

021 

36 

0.91 

022 

32 

0.81 

023 

28 

0.71 

024 

24 

0.61 

025 

20 

0.51 

026 

16 

0.41 

027 

12 

0.30 

028 

10 

0.25 

029 

8 

0.20 

030 

6 

0.15 

031 

4 

0.10 

In  the  generation  of  the  reference  set  and  all  other  images,  one  addi- 
tional operation  was  performed  to  compensate  for  the  specifics  of  sampling, 
filtering,  and  reconstructing  operation.  This  operation  was  a reading  and 
writing  spot  transfer  function  compensation,  in  which  the  Fourier  transform 
of  the  sampled  image  was  divided  by  the  transform  of  the  combination  of  the 
reading  and  writing  spots  used  to  sample  and  reconstruct  the  images.  For  the 
object  images,  this  reading  and  writing  spot  transfer  function  compensation 
effectively  eliminates  the  reading  and  writing  spots  as  sources  of  modulation 
transfer  function  (MTF)  degradation. 
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2.2.2  Generation  of  Short  Time  Exposure  Images 

For  the  generation  of  the  short  time  exposure  transfer  functions,  ran- 
dom wavefronts  were  generated  whose  phase  amplitude  variations  exhibited 
kolmogoroff  statistics,  given  by  a power  spectrum  of  the  form 


W (k)  = 0.023r  “5/3  k"11/3 

(J)  O 


2 2 
[radians  /(cycle/meter)  ] (3) 


where 


V 


2 2 
k + k 
x y 


r = coherence  length  of  the  atmosphere 
o 

These  wavefronts  were  then  Fourier  transformed,  to  obtain  the  optical 
transfer  functions  (OTF)  that  were  subsequently  used  to  generate  the  images. 

An  alternative  description  of  a wavefront  exhibiting  Kolmogoroff  statis- 
tics is  found  by  specifying  the  phase  structure  function,  defined  as 


D(r)  - + r)  - <t>(r')|^) 


(A) 


where 


J 2 2 

r ■ Jl  u + v 

4>(r’)  * phase  of  wavefront  at  point  r' 

r ■ spatial  separation  between  two  points  in  wavefront 

This  wave  structure  function  can  be  written  for  the  case  of  Kolmogoroff  tur- 
bulence as 

/ \5/3 

D(r)  = 6.88(  — j 

where  r describes  the  strength  of  the  atmospheric  turbulence.  Typical  values 
o 

of  the  coherence  length  for  typical  seeing  conditions  are  in  the  range  of  4 
inches  (10  centimeters). 


II 


2. 2. 2.1  Wavefront  Simulation  - There  are  several  methods  available  for  the 
simulation  of  wavefronts  having  phase  perturbations  exhibiting  Kolmogoroff 
statistics:  the  method  of  generating  the  wavefronts  used  in  this  study  in- 
volves describing  the  wavefront  as  a Zernike  polynomial  expansion  given  by 

00 

()>(u,v)  ■ E A^Z^u.v) 
h n“l 


where 

<()(u,v)  * wavefront  perturbation  (waves) 

Z (u.v)  ■ nth  orthonormal  Zernike  polynomial 
n * 

A - coefficient  of  nth  Zernike  polynomial 
n 

Note  that  any  given  wavefront  is  described  by  the  set  { An } and  that  for 
Kolmogoroff  turbulence, 

<\>  ‘ °'°  2 <7) 
and  the  value  of  the  mean  square  quantity,  /Ar  ^),  can  be  computed  from  the 

Kolmogoroff  power  spectrum  and  the  Fourier  transforms  of  the  Zernike  poly- 
nomials. 

Wavefront  simulation  was  performed  by  using  a random  number  generator  to 

select  a set  of  coefficients  {A  },  whose  mean  and  variance  satisfy  the  condi- 

n 

tions  stated  above.  The  advantage  of  this  method  over  the  generation  of  a 
wavefront  by  Fourier  series  representation  lies  in  the  fact  that  each  of  the 
terms  in  the  expansion  is  explicitly  accounted  for;  there  are  no  errors  in 
the  generation  of  the  low  spatial  frequency  terms  (low  order  errors,  such  as 
wavefront  tilt)  resulting  from  finite  sampling  in  the  frequency  domain.  Ap- 
pendix A contains  a detailed  discussion  of  this  wavefront  generation  method. 

Table  3,  below,  lists  expected  values  of  rms  wavefront  error  for  various 

values  of  coherence  length  r , for  two  cases  of  interest  - short  and  long 

o 

time  exposures.  Note  that  the  long  time  exposure  case  includes  the  wavefront 
error  due  to  the  tilt  component,  while  the  short  time  exposure  case  does  not. 
(In  the  case  of  a short  time  exposure,  tilt  causes  only  a displacement  of  the 
object,  which  does  not  degrade  its  quality.) 
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TABLE  3.  WAVEFRONT  ERRORS  FOR  60-INCH  APERTURE 


r 

0 

Coherence 

Length 

(Meters) 

Wavefront 

Error 

Short  Time 

Exposure 

Long  Time  Exposure  j 

Mean  Square 
Error  (Rad2) 

rms  Waves 

Mean  Square 
Error  (Rad2) 

rms  Waves 

1.50 

0.13 

1.00 

0.16 

0.80 

0.39 

3.00 

0.28 

0.40 

1.25 

0.180 

9.61 

0.49 

0.25 

2.74 

0.260 

21.06 

0.73 

0.15 

6.43 

0.40 

49.42 

1.12 

0.10 

12.67 

0.57 

97.39 

1.57 

1 2 

2. 2. 2. 2 Optical  Transfer  Function  Generation  - Both  Fried  and  Hufnagel 
predict  that  the  expected  value  of  the  short  time  exposure  transfer  function 
Tge(k)  for  viewing  through  a turbulent  medium  exhibiting  Kolmogoroff  statis- 
tics is  given  by 


where 

<>  = denotes  an  ensemble  average 

T (k)  =>  diffraction-limited  transfer  function 
o 

X = wavelength 

F * focal  length 

Dq  * aperture  diameter 

r * coherence  diameter 
o 


D.L.  Fried,  "Optical  Resolution  Through  a Randomly  Inhomogeneous  Medium 
for  Very  Long  and  Very  Short  Exposures",  J.  Opt.  Soc.  Amer.,  Vol.  56, 
No.  11,  1372  (1966). 

2 

R. E.  Hufnagel,  "On  the  Mean  Short-Term  Blur",  Woods  Hole  Summer  Study  on 
Restoration  of  Atmospherically  Degraded  Images,  Vol.  2,  Appendix  4 (July 
1966). 


13 


However,  this  MTF  is  not  descriptive  of  the  OTF  of  a single  short  time  expo- 
sure OTF,  since  the  process  of  performing  the  ensemble  average  involves  find- 
ing the  average  magnitude  ^T(k)^  of  a complex  quantity  (the  OTF)  whose  phase 
angle  is  essentially  a random  quantity;  thus,  although  the  MTF  may  have  a sig- 
nificant average  value  independent  of  the  phase,  the  inclusion  of  the  phase 
portion  of  the  OTF  in  the  averaging  process  performed  in  computing  <[jse00^ 
will  cause  the  function  to  drop  to  zero  rapidly. 


Korff3  has  predicted  the  expected  value  of  the  square  of  the  absolute 
value  of  the  short  time  exposure  transfer  function  as 


<lTse«oi2> 


T (k) 
o 

<Do/ro): 


0.435  + 


0.278 


/ D k \1/3 

t 


D k 

for  >>  1 
o 


(9) 


where  k - normalized  spatial  frequency 


Since  this  quantity  is  computed  by  ignoring  the  phase  portion  of  the  OTF, 
it  is  this  quantity  that  should  correspond  to  the  average  modulus  of  simulated 
OTF's  generated  from  random  wavefronts  exhibiting  Kolmogoroff  statistics. 


To  check  the  assumption  stated  above,  a series  of  random  wavefronts  was 
generated  for  an  rQ  of  13  centimeters,  and  used  to  generate  a series  of  opti- 
cal transfer  functions.  Figure  4 illustrates  values  of  the  various  generated 
OTF's,  along  with  an  evaluation  of  the  Korff  expression  for  the  average  short 
time  exposure  MTF.  As  shown,  the  agreement  is  quite  good.  This  model  will 
be  used  later  for  the  purpose  of  deriving  post-processing  filters  based  on 
only  a statistical  estimate  of  the  degrading  transfer  function. 


2.2.3  Image  Generation 

Using  the  software  described  in  Appendix  A,  a series  of  images  containing 
the  simulated  effects  of  various  amounts  of  atmospheric  turbulence  was  gen- 
erated for  the  five  values  of  r listed  in  Table  4.  Some  of  these  transfer 

o 

3D.  Korff,  "Analysis  of  a Method  for  Obtaining  Near-Diffraction  Limited  In- 
formation in  the  Presence  of  Atmospheric  Turbulence",  J.  Opt.  Soc.  Amer., 

Vol.  63,  No.  8,  971  (1973). 
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transfer  functions  are  shown  in  Figure  5.  Paragraph  2.3  contains  a discussion 
of  the  production  of  noisy  images,  in  which  various  amounts  of  both  Poisson 
and  Gaussian  noise  are  added  to  these  atmospherically  degraded  images. 

TABLE  4.  SUMMARY  OF  NOISELESS  SHORT  TIME  EXPOSURE  IMAGES 


Image 

Number 

r 

0 

(Meters) 

Mean  Square 
Wavefront 
Error  (Rad^) 

Strehl 

Ratio 

015 

OO 

0.00 

1.0000 

052 

0.80 

0.39 

0.6900 

048 

0.40 

1.25 

0.3400 

044 

0.25 

2.74 

0.1300 

040 

0.15 

6.43 

0.0320 

036 

0.10 

12.67 

0.0044 

2.2.4  Generation  of  Long  Time  Exposure  Images 

For  the  generation  of  long  time  exposure  images,  the  assumption  will  be 
made  that  the  images  will  be  obtained  through  a pre-compensation  system  in 
which  a significant  amount  of  wavefront  correction  is  performed  before  the 
image  is  recorded,  with  the  emphasis  on  the  removal  of  the  lew  spatial  fre- 
quency wavefront  errors. 

Hulnagel  and  others  have  shown  that  the  expected  value  of  the  long  time 
exposure  MTF,  for  no  wavefront  correction,  is  given  by: 


<^T£e Ck^>  = To(k)e'[C(o)  ' C(XFk)] 


where  C(AFk)  * C(r)  and  is  the  two-dimensional  autocorrelation  function  of 

2 

the  wavefront  deformations,  in  units  of  radians  . The  phase  structure  func- 
tion, D^(r),  is  given  as 

Vr)  = ^l$<r  + r'>  ~ <Kr)|2)>  (11) 
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Expanding  Equation  (11)  yields 


D^r)  = 2[C(0)  - C(r) ] (12) 

Equation  (10)  becomes,  upon  substitution  of  Equations  (5)  and  (12) 

-3.44  imiyn 

(Tlle(k))  =To(k)e  ^ <13> 

For  the  case  where  a significant  amount  of  wavefront  correction  has  been  per- 
formed such  that  the  low  spatial  frequency  components  of  the  wavefront  have 
been  reduced,  the  autocorrelation  function  of  the  wavefront  will  have  a lower 
value  at  AFk  = 0,  and  will  drop  faster,  resulting  in  a phase  structure  func- 
tion that  levels  off  at  large  values  of  AFk  (see  Figure  6).  For  this  case, 
the  long  time  transfer  function  becomes  approximately 

A.  (kj)  = T (k)c‘1/2D(XFk)  * T (k)e"C(o)  (14) 

\ Joe  / o o 

where  C(o)  is  the  mean  square  uncorrected  wavefront  error. 

Figure  7 illustrates  a family  of  four  long  time  exposure  transfer  func- 
tions, computed  from  Equation  (14)  for  which  an  approximate  shape  was  assumed 
for  the  spatial  frequency  region  near  the  origin;  Table  5 lists  the  noise- 
less images  generated  from  these  transfer  functions;  and  Figure  8 illustrates 
three  of  the  noiseless  long  time  exposure  images. 


TABLE  5.  SUMMARY  OF  LONG  TIME  EXPOSURE  IMAGES 


Image 

Number 

Mean  Square 
Wavefront  Error  (Rad  ) 

Strehl 

Ratio 

089 

0.39 

0.838 

090 

1.20 

0.391 

091 

2.74 

0.109 

092 

4.36 

0.031 

No  Wavefront  Correction 
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Figure  6.  Wavefront  Autocorrelation  Functions 
and  Phase  Structure  Functions 


Note  that-  a direct  comparison  of  short  and  long  time  exposure  images  reveals 
that  the  short  time  exposure  images  differ  from  the  long  time  exposure  images 
by  having  much  more  fine  structure  visible,  because  of  the  generally  higher 
Modulus  values  and  the  presence  of  substantial  phase  in  the  OTF. 


2.3  NOISE  SIMULATION 

The  second  important  image  quality  degrading  feature  of  the  Compensated 
Imaging  System  to  be  simulated  is  image  noise.  There  are  two  basic  types  of 
image  noise  that  the  Compensated  Imaging  System  may  be  subject  to  - signal- 
dependent  (shot)  noise  arising  from  the  statistics  of  photon  arrival  rates, 
and  additive  noise.  Both  of  these  noise  mechanisms  have  been  simulated  in 
the  study,  as  described  below. 


2.3.1  Poisson  (Shot)  Noise  Simulation 

In  the  image  recording  process,  an  image  sampling  operation  is  per- 
formed, resulting  in  a sampled  image  consisting  of  a series  of  digital  num- 
bers representing  the  number  of  photoelectrons  collected  from  each  image 
picture  element  (pixel)  during  an  exposure  (photoelectron  integration  time). 
The  value  of  that  sampled  image  data  point  may  be  written  in  terms  of  the 
reflectivity  of  the  object  and  a photometric  constant  that  lumps  together 
all  of  the  optical  efficiency,  collecting  aperture,  and  object  irradiance 
factors. 


where 


U (x)  *=  A r (x)  (photoelectrons)  (15) 

s s 


U (x)  * value  of  sampled  image  (photoelectrons) 
s 

A ■ photoelectrons  per  pixel  per  unit  object  reflectivity 

l (x)  ■ object  reflectivity 
s 

The  rms  shot  noise  On(x)  resulting  from  the  signal  described  in  Equation 
(15)  can  be  written  as 


a (x)  * 4/U  (x)  n(x)  (photoelectrons)  (16) 
n ■ s 


(17a) 


where  n(x)  is  a random  variable  exhibiting  Poisson  statistics,  for  which 
^n(x)^>  - 0.0 

(in(x)l2)  -1.0  (17b) 

We  may  now  express  the  rms  noise  in  equivalent  reflectivity  units  by  using 
the  constant  A 

O (x)  / r (x) 

a (x)  _ _S - %l-~ — n(x)  (reflectivity)  (18) 

r A " A 

Taking  an  average  over  the  encire  satellite  yields 

(reflectivity)  (19) 

For  the  simulation  of  Poisson  noise,  Equation  (18)  was  used  in  conjunc- 
tion with  an  algorithm  that  permits  rapid  generation  of  values  of  the  Poisson 
random  variable  defined  by  n(x).*  Values  of  A used  ranged  from  16  photoelec- 
trons per  pixel  per  unit  reflectivity  to  2000  photoelectrons  per  unit  reflec- 
tivity. This  range  of  values  of  A resulted  in  a range  of  average  rms  reflec- 
tivity values  from  0.0112  reflectivity  unit  to  0.125  reflectivity  unit,  based 
on  an  average  reflectivity  of  25  percent.  * 

2.3.2  Poisson  and  Gaussian  Noise  Image  Generation 

Five  secs  of  noisy  images  were  generated  for  the  study,  as  summarised  in 
Figure  1.  The  first  two  sets  (B  and  H),  consisting  of  noise  images  obtained 
through  diffraction-limited  apertures  of  60-,  48-,  32-,  and  16-inch  diameters, 
were  generated  to  provide  images  against  which  various  image  quality  models 
could  be  measured.  The  other  three  sets  were  generated  to  provide  the  type 
of  noisy  Imagery  that  a Compensated  Imaging  System  might  be  expected  to  pro- 
duce, for  the  case  of  pre-compensation  only  (the  case  of  post-compensation 
is  discussed  in  paragraph  2.4).  See  Figure  9 for  example. 

Figures  10  through  14  illustrate  the  selection  of  the  0TF  and  noise 
parameters  for  each  of  the  simulations;  the  three-digit  numbers  in  each  of 
the  rectangles  refers  to  the  frame  number  of  that  specific  simulation. 

4R.E.  Hufnagel  and  E.L.  Kerr,  "A  Simple  Algorithm  for  Fast  Real-Time 
Generation  of  Pseudorandom  Poisson  Integers  with  Rapidly  Varying  Means  , 
Proc.  IEEE,  Vol.  57,  No.  11,  2088  (1969) 
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a - 80  a « 16 

- 0.056  or  - 0.125 


Figure  9.  Noise  Simulation  Examples 
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Figure  13.  Set  I.  Short  Time  Exposure  Atmospheric 
OTF  and  Gaussian  Noise 


28 


*• 


Appendix  B contains  nine  sample  images  from  each  of  the  five  sets  summarized 
in  Figures  10  through  14. 

2.4  POST-DETECTION  COMPENSATION  SIMULATION 

One  of  the  basic  features  or  the  Compensated  Imaging  System  concept  is 
the  use  of  post-detection  image  processing  to  perform  additional  image  res- 
toration, after  the  pre-detection  compensation  adaptive  optics  has  done  as 
much  as  possible  to  correct  the  wavefront  to  that  of  a diffraction-limited 
system. 

If  image  restoration  is  accepted  as  the  prime  goal  of  the  system,  then 
the  post-detection  compensation  operation  has  as  its  goal  the  restoration 
of  the  system  transfer  function  to  that  of  a diffraction-limited  aperture 
having  the  same  aperture  diameter  as  the  system  itself.  Because  of  the 
presence  of  noise  in  the  images  being  restored,  the  restoration  goal  may  not 
yield  the  most  desirable  post-detection  compensated  image  in  all  cases;  how- 
ever, this  goal  of  restoring  the  transfer  function  is  the  most  straight- 
forward, and  is  likely  to  provide  the  basis  for  most  of  the  operations  to  be 
performed  on  imagery  having  good  signal- to-noise  ratios  in  any  operational 
system.  In  addition,  for  cases  where  the  system  is  not  noise  limited,  the 
results  of  such  post-detection  transfer  function  compensation  will  appear  to 
be  equal  to  images  obtained  from  a diffraction-limited  system,  which  is  the 
basic  goal  of  the  system  concept. 

There  are  two  considerations  to  be  made  in  the  simulation  of  the  post- 
detection compensation  function;  first,  what  is  the  effect  of  imperfect  know- 
ledge of  the  degrading  transfer  function  and  second,  what  is  the  effect  on 
the  subjective  quality  of  the  changed  spectrum  of  the  image  noise  due  to  the 
effects  of  application  of  post-detection  compensation,  compared  to  the  ap- 
plication of  the  same  general  level  of  compensation  applied  before  image 
detection. 

2.4.1  Exact  Transfer  Function  Restoration 

If  the  exact  optical  transfer  function  through  which  an  image  is  obtained 
is  known,  the  post-detection  compensation  operation  can  take  the  form  of  a 
direct  inverse  frequency  domain  filter,  as  described  below,  where  the  filter 


JO 


; 


■ 


F (k)  is  derived  from  the  known  degrading  transfer  function  and 

the  diffraction-limited  transfer  function  T^k) . (The  values  of  E^^  and 
3 serve  to  avoid  singularities  in  the  T function.) 

c B 


F (k) 
e 


0 

TJ0(k) 


Tkn(k) 

TJ0(k) 


Wk> 

^max 

> T. 
kn 

E 

max 

\„<k) 

T 

< ^ 

(k)  > 
(k) 


(20) 


max 


For  this  simple  type  of  restoration  operation,  the  parameter  Emax  can 
be  varied  to  produce  different  post-processing  Strehl  ratios.  For  the  images 
contained  in  sets  D,  K,  and  G,  which  represent  the  application  of  this  type 
of  restoration  filter  to  the  short  and  long  term  exposure  image  sets  C,  I, 
and  F.  the  value  of  E for  each  of  the  post-detection  filters  was  set  to 

produce  a Strehl  ratio  for  the  post-processed  images  of  0.8.  Figure  15  il- 
lustrates the  three  post-proct-ssed  transfer  functions  and  corresponding 
images,  for  the  post-detection  compensation  of  the  short  time  exposure  im- 
ages. The  examples  shown  are  the  Poisson  noise  images.  Figures  16  and  17 
illustrate  the  combination  of  parameters  for  which  th  exact  restoration 
filters  illustrated  in  Figure  15  have  been  applied  to  the  Poisson  and 
Gaussian  noise  degraded  short  time  exposure  images;  Figure  18  illustrates 
the  parameters  for  the  analogous  operation  applied  to  the  long  time  expo- 
sure images,  for  the  case  of  Poisson  noise  degradation. 

2.4.2  Approximate  Transfer  Function  Restoration 

In  addition  to  the  exact  restoration  filter  discussed  above,  a second 
type  of  restoration  filter  can  be  used,  one  based  on  a statistical  know- 
ledge of  the  transfer  function,  rather  than  on  exact  knowledge.  The  Korff 
formulation  for  the  average  modulus  of  a short  time  exposure  transfer  func- 
tion was  discussed  above;  if  that  formulation  is  used  as  the  basis  for  de- 
riving an  inverse  filter,  as  described  below,  the  result  of  applying  that 
filter  to  a short  time  exposure  image  will  be  a post-processing  transfer 
function  in  which  some  residual  modulus  and  phase  error  will  remain 
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Figure  IS.  Short  Tine  Exposure  OTF  end  Post-Processing  (Exact) 
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Figure  17.  Set  K.  Short  Time  Exposure  Atmospheric  OTF 
and  Gaussian  Noise  and  Post-Processing  to 
Strehl  Ratio  ■ 0.80 
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Figure  18.  Set  G.  Long  Time  Exposure  JITF  and  Poisson 
Noise  and  Post-Processing  Restoration  to 
Strehl  Ratio  ■ 0.8 
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See  Figure  19  for  three  examples ; Figure  20  summarizes  the  images  generated . 

As  shown  in  Table  6,  which  lists  the  pertinent  parameters  for  all  of  the  post- 
processed  images,  the  application  of  the  approximate  filter  results  in  Strehl 
ratios  that  fall  considerably  short  of  the  Strehl  ratio  for  the  post-processed 
images  processed  using  exact  knowledge  of  the  degrading  transfer  function. 


TABLE  6.  SUMMARY  OF  POST-PROCESSED  TRANSFER  FUNCTION 
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Figure  19.  Short  Time  Exposure  OTF  and  Post-Processing 
(Approximate) 
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Figure  20.  Set  J.  Short  Time  Exposure  Atmospheric  OTF  and  Poisson 
Noise  and  Post-Processing  Using  Korff  Model 


SECTION  III 


VIEWING  AND  INTERPRETATION  EXPERIMENTS 


In  Section  I,  the  main  purpose  of  the  study  described  in  this  report 
was  stated  to  be  a derivation  of  a mathematical  model  relating  the  subjec- 
tive quality  of  satellite  imagery,  as  determined  by  human  observers,  to 
various  quantitative  measures  of  system  transfer  function  and  noise  charac- 
teristics. In  this  section,  the  determination  of  the  subjective  image  qual- 
ity of  the  Images  by  human  observers  will  be  discussed. 

3.1  VIEWING  AND  INTERPRETATION  EXPERIMENT  DESIGN 

The  most  direct  method  of  determining  the  relationship  between  the  sub- 
jective image  quality  and  the  degradations  present  in  the  satellite  Imagery 
is  to  perform  viewing  and  interpretation  experiments.  In  these  experiments, 
a group  of  human  observers  is  asked  to  rank  each  of  the  degraded  images 
.against  a reference  set  of  standard  Images  having  known  resolution  charac- 
teristics. If  the  observer  decides  that  a degraded  image  is  subjectively 
equal  to  one  of  the  standard  set,  then  the  degraded  image  is  said  to  have 
an  image  quality  equal  to  t&e  quality  of  the  standard  image  to  which  it  was 
equated  (it  is  also  permissible  for  the  observer  to  rate  an  image  between 
two  members  of  the  reference  set) . Previous  studies  have  shown  a high  cor- 
relation between  the  subjective  quality  judgements  of  a group  of  observers, 
and  the  actual  performance  of  professional  photo-interpreters  measured  by 
the  amount  of  quantitative  information  correctly  recovered  in  a controlled 
experiment 

For  the  experiments  performed  in  this  study,  the  resolution  of  the  image, 
as  determined  by  the  average  relative  ranking  of  that  image  against  the  ref- 
erence set,  is  expressed  in  terms  of  equivalent  aperture  diameter.  Thus,  an 
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60-inch-aperture  system  may  well  have  the  resolution,  or  equivalent  quality, 
of  a diffraction-limited  image  obtained  from  a 4-inch-aperture  system.  As 
will  be  seen  in  Section  IV,  these  units  for  expressing  image  quality  are 
compatible  with  several  image  quality  measures  of  interest. 

All  of  the  simulations  performed  for  atmospherically  degraded  transfer 
functions  were  generated  for  a 60-inch-aperture  system;  only  sets  A,  B,  and  H 
contained  images  generated  for  apertures  smaller  than  60  inches. 

i 

3.2  DESCRIPTION  OF  RESULTS  OF  THE  V AND  I EXPERIMENTS 

The  106  images  produced  for  the  study  were  separated  into  seven  separate 
groups,  and  each  group  of  images  was  viewed  by  an  average  of  13  observers. 

Each  observer  was  asked  to  rank  each  of  the  images  against  the  reference  set 
of  15  images  (see  Figure  21) . 

For  each  image,  the  geometrical  mean  of  the  quality  estimates  was  de- 
termined, along  with  the  standard  deviation  factor  for  the  set  of  estimates. 

The  distribution  of  the  estimates  for  any  single  image  was  found  to  be  ap- 
proximately a log  normally  distributed  Gaussian  random  variable;  therefore, 

I ! j 

the  geometrical  mean  should  give  the  truest  estimate  of  the  average  per- 
ceived quality  for  each  of  the  images.  Figure  22  illustrates  typical  cumula- 
tive probability  plots  of  the  individual  quality  estimates  for  three  of  the 
inages  viewed. 

Tables  7 through  9 list  the  resultant  estimates  of  subjective  image 
quality  determined  for  each  of  the  106  images.  Table  7 lists  the  perceived 
quality  for  each  of  the  10  short  and  long  term  exposure  noiseless  images; 

Table  8 lists  the  results  for  all  of  the  Poisson  noise  images;  and  Table  9 
lists  the  results  for  all  of  the  Gaussian  noise  images.  For  the  tables  shown, 
the  column  headings  have  the  following  meaning: 

IMAGE  - image  frame  number 

DIA  - aperture  diameter  in  inches 
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Figure  21.  Viewing  and  Interpretation  Experiment 


Figure  22.  Cumulative  Probability  Distributions  of  Quality  Estimates 
for  Three  Typical  Images 
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TABLE  8.  POISSON  NOISE  IMAGE  V AND  I RESULTS 
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TABLE  8.  POISSON  NOISE  IMAGE  V AND  I RESULTS  (CONTINUED) 
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TABLE  9.  GAUSSIAN  NOISE  IMAGE  V AND  I RESULTS 
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- the  factor  by  which  the  standard  deviation  of  ad- 
ditive Gaussian  noise  has  been  increased  by  the 
application  of  po*>l  -detection  compensation  fre- 
quency domain  filter 

the  Strehl  ratio  for  the  particular  OTF,  normal- 
ized to  the  diffraction-limited  Strehl  ratio  for 
that  aperture;  for  the  case  of  post-detection 
processed  Images , the  Strehl  ratio  is  the  post- 
processed  Strehl  ratio 

the  value  of  the  noise  generation  parameter 
(photoelectrons  per  pixel  per  unit  reflectivity) 

- the  value  of  the  rms  noise  in  reflectivity  units, 
before  any  post-detection  processing 

- the  value  of  the  average  quality  of  the  image,  as 
determined  by  the  observers  (units  are  inches  of 
equivalent  aperture) 

standard  deviation  factor  for  the  set  of  quality 
estimates 

IV,  the  fitting  of  this  raw  data  to  various  image  quality 


measures  will  be  discussed. 


SECTION  IV 


ANALYSIS  OF  RESULTS 


4.1  IMAGE  QUALITY  CRITERIA 

The.  e are  several  classes  of  image  quality  criteria  that  are  available, 
as  shown  in  Table  10;  the  entries  in  Table  10  are  listed  in  terms  of  the 
information  concerning  the  image  formation  process  that  is  used  in  the  evalu- 
ation of  each  image  quality  criteria.  It  is  to  be  expected  that  the  measures 
that  use  the  most  information  will  be  the  ones  that  describe  system  perform- 
ance. 

4.1.1  Class  A:  Strehl  Ratio,  Edge  Sharpness,  and  Normalized  Relative  Energy 

These  measures  are  all  characterized  by  the  fact  that  each  is  a function 
only  of  the  system  transfer  function.  The  Strehl  ratio  describes  the  value 
of  the  point  spread  function  at  the  centroid,  and  is  given  by 

fT(k)dk 

s - * p(o,o^  t22) 

/Idl(k)dk 


T(k)  = system  optical  transfer  function 

T (k)  = diffraction-limited  transfer  function  for  same  size 

Q X< 

aperture 

p(0,0)  = value  of  point  spread  function  at  the  centroid 

The  unnormalized  Strehl  ratio  has  the  units  of  spatial  frequency  squared; 
thus,  the  square  root  of  the  Strehl  ratio,  unnormalized,  has  the  units  of 
spatial  frequency,  and  represents  an  effective  spatial  frequency  bandwidth 
for  the  system  whose  transfer  function  is  given  by  T(k).  In  addition,  the 
square  root  of  the  unnormalized  Strehl  ratio  is  proportional  tc  the  diameter 
of  the  aperture  of  the  system;  therefore,  it  lends  itself  dimensionally,  at 
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TABLE  10.  IMAGE  QUALITY  MEASURES 


Information  Used 

Type 

Image  Quality  Measure 

Transfer 

Function 

Noise 

Object 

Display 

A 

Strehl  Ratio 

Edge  Sharpness 

Normalized  Relative  Energy 

X 

B 

Image  Sharpness 
Image-Object  Correlation 
Image  Fidelity 

X 

X 

C 

Summary  Measure 

X 

X 

X 

least,  to  use  as  a measure  of  image  quality  having  the  units  of  equivalent 
aperture  diameter.  However,  the  use  of  the  Strehl  ratio  as  a measure  of 
image  quality  sometimes  breaks  down  for  comparison  of  images  obtained  through 
the  transfer  functions  of  widely  differing  shapes. 

The  edge  response  of  the  system  describes  the  slope  of  the  edge  response 
of  the  system,  and  is  a measure  of  how  sharply  edges  are  imaged.  For  systems 
that  contain  edges,  this  measure  may  be  more  appropriate  than  the  Strehl 
ratio,  which  describes  the  imaging  properties  for  point  objects.  The  edge 
response  is  given  by 


E 


£ “W1 


X“0 


- j W> dk* 


(23) 


where 

E “ slope  of  edge  response 
i,(x)  - line  spread  function 


The  normalized  relative  energy  (NRE)  criterion  describes  the  fraction 
of  the  total  energy  contained  in  the  system  spread  function  that  is  contained 
within  a circle  of  some  diameter  centered  on  the  centroid  of  the  point  spread 
function.  This  criterion,  similar  to  the  Strehl  ratio,  is  most  suitable  for 
systems  in  which  detectors  are  used  to  make  photometric  measurements.  The 
NRE  is  given  by 


J (2TTRk) 

NRE  - J p(r)dr  - /l(k)  ^rrT  dk  <24> 

A-TTR2 


4.1.2  Class  B:  Image  Sharpness,  Image-Object  Correlation,  Mean-Square  Error 

This  class  of  image  quality  measures,  representing  two-dimensional  gen- 
eralizations of  statistical  measures  used  for  evaluating  time  filters,  was 
introduced  by  Linfoot^,  and  is  discussed  by  O'Neill2. 


The  image  sharpness  criterion,  which  describes  the  relative  structural 
content  of  an  image,  relative  to  the  structural  content  of  the  object,  is 
given  by 


t2(r)  J* | T (k ) ]2  $(k)dk 

o2(r)  J^OOdk 


(25) 


where 


i(r)  * image  in  two  dimensions 
o(r)  = object 

<t>(k)  = power  spectrum  of  object 

The  image-object  correlation  measure  is  given  below,  and  measures  the 
alignment  of  the  object  and  the  image: 


^E.H.  Linfoot,  J.  Opt.  Soc.  Amer.  , 46 , 740  (1956). 

2E.  O'Neill,  Introduction  to  Statistical  Optics,  Addison  Wesley,  1963,  p.  106. 
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m 


ioc 


|*T(k)*(k)dk 
i(r)o(r)__  . J 

o2(r)  J *tV.)dk 


(26) 


The  meSEJBHJIS-HISS.  b«~»  the  *n'1  th*  0bleCt  “ 8iV“  by 

f^T-T(J)]2  _ f[l-2I(k)4  |T(k)|21«(k)dk 

D * " , / 1 \ jl. 


(27) 


1*.SC 


o2(r) 


T t(k)dk 


Like  the  Strehi  retie,  the  edge  reeponee,  end  the  NRE,  theee  three  ■“»“  ' 

ere  not  in  theeselve,  euitebl.  for  predicting  the  quelity  of  i-ge.  for 
tens  where  there  (e  eignlfioent  imege  noiee. 

4 1.3  dess  C!  Sumnery  Meeeure  of  Imege  Quelity 

Based  on  eppro.i-.tely  ten  yeer.  of  inveetlgetlon.  into  the  netor.  of 
i-ege  quelity.  Dr.  R.E.  Hufnegel  of  P«rkln-gl»er  he.  for-uUted  • «—» 
-sure  of  1— ege  quelity  thet  teke.  into  eccount  the  ey.te.  tr.ne««r  ’ 

the  1— ege  noUe  present  in  the  systen,  the  transfer  function  of  the  di  P y 
medium  used  to  present  the  imege  to  the  viewer  (or  po.t-proc.ee  the  i-eg  ). 
and  the  transfer  function  of  the  eye  of  the  observer. 


opt 


(28) 
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(k  »k  ) measured  at  the  display  plane 
x y 

dk  dk 
x y 

system  OTF  up  to  the  point  where  noise  is  introduced 

display  OTF  including  image  restoration  attempts  (after 
noise  introduction)  and  display  spot  shape  effects 

Nyquist  frequency  bounds,  when  the  image  is  spatially  quan- 
tized. In  the  absence  of  sampling 

human  visual  system  MTF  for  25-cm  viewing  distance,  with 
k in  units  of  cy/mm 

magnification  of  viewing  ( *=  — — — -cm \ 

\ viewing  distance  for  unaided  eye) 

noise  Wiener  spectral  density  for  additive  Gaussian  noise. 

If  the  noise  is  not  additive  Gaussian,  the  closest  equiva- 
lent is  used  (e.g.,  Poisson  noise  is  evaluated  at  mean  re- 
flectivity). Note  that  the  noise  is  to  be  expressed  in  units 
of  (object  reflectivity)2/ (cy/mm)2.  Expressing  the  noise  as 
an  apparent  fluctuation  of  object  reflectivity  circumvents 
problems  of  illumination  levels,  display  brightness,  camera 
nonlinearities,  etc.  The  eye-brain  system  automatically 
compensates  for  all  these  factors.  If  spatial  sampling  is 
used,  the  folded  over  noise  spectrum  is  to  be  used. 

noise  scaling  factor  (to  be  discussed) 

pre-noise-introduction,  modulation  ratio  factor  to  account  for 
change  of  modulation  due  to  haze  or  sensor  gamma.  No  change 
implies  M = 1 . 

post-noise  introduction  modulation  factor  including  display. 

It  is  assumed  that  neither  M nor  is  big  enough  to  cause 

significant  loss  of  information  by  nonlinear  saturation  ef- 
fects. 

overall  image  quality  in  units  of  spatial  frequency 

dimensionless  scale  factor  (if  this  factor  is  the  aperture 
diameter,  Q has  the  units  of  equivalent  aperture  diameter). 
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The  right-hand  portion  of  the  denominator  is  the  mean  squared  noise  as  seen 
with  the  saae  magnification.  The  factor  a is  a noise  scaling  term  with  the 
units  (object  reflectivity)  ^ and  should  be  associated  with  the  mean  squared 
modulation  in  the  objec  scene.  For  aerial  photographs  of  ground  terrain, 
a ■ 1600  has  been  found  to  give  a best  fit  to  experiments. 

4.2  EXPERIMENTAL  DATA  FIT  TO  SIMPLE  SUMMARY  MEASURE  MODEL 

The  summary  measure  of  image  quality  model  contains  a factor  ct,  which 
has  not  been  determined  for  imagery  of  the  type  being  simulated  on  this  study  - 
i.e.,  high-contrast  Images  of  small  objects,  as  opposed  to  extended  images  of 
large  low-contrast  objects  (aerial  scenes)  for  which  there  is  10  years  of 
data.  Thus,  it  was  of  Interest  to  take  a simple  version  of  the  summary  meas- 
ure model,  without  the  eye  MTF,  and  determine  the  value  of  o that  best  fits 
the  experimental  data  from  the  V and  I experiments. 

Once  that  value  of  ot  has  been  determined , then  the  model  can  be  used  to 
estimate  the  quality  of  each  of  the  Images  viewed  in  the  study;  the  rms  error 
between  the  estimates  made  on  the  basis  of  the  a value  and  the  actual  observed 
values  will  indicate  the  ability  of  the  model  to  handle  the  substantial  varia- 
tion conditions  for  which  the  simulations  were  generated.  Once  this  step  is 
complete,  then  the  full  model  can  be  Investigated. 

Equation  (28)  describes  the  model  for  the  case  where  the  eye  MTF  has 
been  set  to  unity;  as  shown,  the  measure  reduces  to  a simple  combination  of 


J*  T(k)T,  (k)dk 


.1/2 


^l  + otNo  J*|Td(k)|2dk  <1  + a0rNI  ) 


074 


where 


(29) 


D “ aperture  diameter 
o 

2 2 
a - noise  variance  (reflectivity  ) 
r 

Nj  - additive  Gaussian  standard  deviation  increase  factor 
(j  1 when  post-processing  is  performed) 
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4.2.1  Fit  of  Noiseless  Images  to  Simple  Model 


Note  that  the  simple  model  given  in  Equation  (29)  reduces  t<>  the  square 
root  of  the  Strehl  ratio  for  the  no  noise  case;  Figure  23  illustrates  the 
fit  of  experimentally  determined  quality  estimates  to  the  predicted  quality 
determined  from  Equation  (29),  evaluated  for  no  noise.  Each  image  viewed  is 
plotted  with  its  abscissa  coordinate  being  the  predicted  quality,  and  the  or- 
dinate being  the  observed  experimentally  determined  quality,  thus  a perfect 
fit  of  the  model  to  the  data  would  result  in  all  points  lying  on  a straight 
H-ne*  As  shown,  the  data  fits  the  straight  line  quite  well  — the  rms  error 
between  the  observed  quality  and  the  predicted  quality  is  2.9  inches  of  equiv- 
alent aperture,  for  all  10  noiseless  Images  generated  for  short  and  long  time 
exposure  OTF's. 

4.2.2  Fit  of  Poisson  Noise  Images  to  Simple  Model 

All  of  the  Poisson  noise  images  (67  in  number)  were  taken  as  a group,  and 
the  a value  that  minimized  the  mean  square  error  between  the  predicted  values 
and  the  observed  values  was  determined  to  be  a value  of  1800.  This  value  of 
a was  then  used  to  derive  predicted  values  for  all  of  the  Poisson  noise  im- 
ages, with  the  resultant  fit  displayed  in  Figure  24.  For  the  value  of  1800, 
the  root  mean  square  error  between  the  predicted  quality  and  the  observed 
quality  is  3.5  inches  of  equivalent  aperture. 

Appendix  C contains  similar  plots  of  the  Poisson  noise  image  results,  but 
for  each  of  the  basic  experiments  plotted  separately.  Close  examination  of 
the  individual  plots  reveals  systematic  deviations  from  the  perfect  straight- 
line  match  between  predicted  and  observed;  some  of  these  deviations,  although 
not  large,  will  be  further  examined  in  a later  section. 

4.2.3  Fit  of  Gaussian  Noise  Images  to  Simple  Model 

All  of  the  Gaussian  noise  images  (30  in  number)  were  taken  as  a group, 
and  the  a value  that  minimized  the  mean  square  error  between  the  observed 
quality  and  the  predicted  quality  was  determined  to  be  a value  of  3400.  Fig- 
ure 25  illustrates  the  degree  of  fit  to  be  about  the  same  as  it  was  for  the 
Poisson  noise  images;  for  the  Gaussian  images,  the  rms  error  between  the 
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Figure  23.  Noiseless  Images  Fit  to  Simple  Model 
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Figure  24.  Poisson  Noise  Images  Fit  to  Simple  Model 
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Figure  25.  Gaussian  Image  Fit  to  Single  Model 
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predicted  and  observed  is  3.3  inches  of  equivalent  aperture  diameter.  Appen- 
dix C illustrates  the  plots  of  the  individual  experiments  performed  for  the 
Gaussian  noise  simulations. 

4.2.4  Conclusions  Concerning  Fit  to  Simple  Model 

As  a result  of  the  application  of  the  simple  model  to  the  total  number  of 
106  images  generated  on  the  study,  the  simple  model  was  found  to  predict  quite 
well  the  observed  quality  of  the  images,  for  a wide  range  of  noise  and  OTF 
conditions.  Therefore,  it  is  concluded  that  the  simple  model  will  serve  well 
as  a tool  for  estimating  the  quality  of  satellite  images,  for  the  type  of 
generic  degradations  simulated  in  this  study. 

It  should  also  be  mentioned  that  the  model  has  a larger  percentage  error 
in  the  region  of  poor  image  quality,  and  performs  on  a percentage  basis  far 
better  for  the  case  of  good  image  quality,  which  is  precisely  the  region  where 
the  tool  will  be  most  needed  - i.e.,  for  tradeoffs  concerning  the  performance 
of  a compensated  imaging  system  that  is  producing  images  not  too  far  from  the 
quality  of  diffraction-limited  images. 

4.3  EXPERIMENTAL  FIT  TO  COMPLETE  SUMMARY  MEASURE  MODEL 

Once  the  examination  of  the  ability  of  the  simplified  model  to  fit  the 
data  was  completed,  the  evaluation  of  the  full  summary  measure  model  includ- 
ing the  eye  MTF  optimization  was  performed.  For  this  case.  Equation  (28) 
becomes 


where  N = noise  spectral  density  for  equivalent  Gaussian  noise.  Note  that 
o 

the  scale  factor  outside  the  square  root  normalizes  the  resultant  values  of 
image  quality  such  that  for  no  noise  Q has  a value  of  60  inches,  for  a 
diffraction-limited  transfer  function  T(k)  of  a 60-inch  aperture  (the  effect 
of  optimizing  the  viewing  magnification  m for  viewing  a diffraction-limited 
transfer  function  is  contained  in  the  0.7845  factor). 
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4.3.1  Fit  of  Noiseless  Images  to  Complete  Model 

For  each  of  the  10  noiseless  images.  Equation  (30)  was  used  to  determine 
the  optimum  (maximum)  value  of  Q.  The  10  resultant  predicted  quality  values 
can  be  compared  to  thij  observed  values  as  determined  by  the  V and  I experi- 
ments, as  shown  in  Figure  26.  The  rms  error  for  the  10  noiseless  imagas  is 
2.0  Inches. 


4.3.2  Fit  of  Noiseless  Images  to  Various  Models 

Table  11  illustrates  a comparison  between  four  measures  of  image  quality, 
again  for  the  noiseless  images.  Shown  are  estimated  quality  values  based  on 
the  full  summary  measure,  the  Strehl  ratio  (equivalent  to  the  simple  sunmary 
measure  model),  and  two  other  figures  of  merit,  as  defined  below). 


4. 3.2.1  Relative  Image  Sharpness  - In  paragraph  4.1,  the  image  sharpness 
criterion  was  given  as 


T l2(r)  J*|T(k)|2*(k)dk 

is  * _ " — : 

o2(r)  I *«lk 

where 

4>(k)  - object  power  spectrum 


(31) 


If  the  object  spectrum  is  assumed  to  be  flat,  and  T is  normalized  with 
respect  to  its  value  for  a diffraction-limited  transfer  function  Tdi(k),  a 
measure  of  relative  image  sharpness  t,  results 

riT(k)i2dk 

t y—  (32) 

,rlTdi<k>!  dk 

This  expression  emphasizes  the  sharpness  of  the  image,  as  measured  by  the  mean 
square  critericn;  there  is  no  explicit  requirement  that  the  image  look  like 
the  object. 
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Noiseless  Images  Fit  to  Full  Model 
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4. 3. 2. 2 Fidelity  Defect  - The  mean  square  error  between  the  image  and  the 
object  was  described  in  paragraph  4.1  as 


D 

mse 


|o(r)  - i(r) |2 
o2(r) 


(33) 


Linfoot  suggests  that  a fidelity  criterion  F,  given  by 


F 

c 


1 - D 

mse 


(34) 


be  used  as  a measure  of  system  quality  (this  measure  behaves  the  same  way  as 
do  the  others  discussed  in  this  report;  i.e.,  it  has  a value  of  unity  for  no 
defect).  If  O'Neill's  work  is  referenced  and  a flat  object  spectrum  is  as- 
sumed, the  equation  can  be  written  as 


where 


f = ar_nt 

l-n 


q = Strehl  ratio  (normalized) 
t = relative  image  sharpness 


and 


n - 


jiTdI(k)i2 

zJlT^OOldk 


(35) 


(36) 


As  shown  in  Table  11,  this  fidelity  criterion  performs  reasonably  well  as 
a measure  of  image  quality;  however,  the  results  of  the  full  summary  measure 
yield  the  best  fit  to  the  results  of  the  V and  I experiments  for  noiseless 
imagery. 


4.3.3  Fit  of  Gaussian  Noise  Images  to  Full  Summary  Measure  Model 

In  c-der  to  fit  the  Gaussian  noise  images  to  the  full  summary  model,  a 
double  optimization  procedure  had  to  be  performed.  The  objective  of  the  pro- 
cedure was  to  determine  the  value  of  a that  minimized  the  mean  squared  error 
between  the  predicted  values  and  the  observed  values  for  all  of  the  Gaussian 
noise  images,  for  the  case  where  the  eye  magnification  was  optimized  for 
each  image  and  value  of  a individually.  The  results  of  this  procedure,  which 
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has  been  implemented  in  a software  package,  is  shown  in  Figure  27;  the  rms 
error  between  the  observed  and  predicted  values  of  image  quality  is  3.5 
inches  rms,  which  is  slightly  higher  than  the  3.3-inch  rms  error  that  re- 
sulted from  the  use  of  the  eye  MTF  Independent  model  described  in  paragraph 
4.2.  The  value  of  a that  produced  the  fit  shown  in  Figure  27  is  2000. 

4.3.4  Fit  of  Poisson  Noise  Images  to  Full  Summary  Measure  Model 

Figure  28  illustrates  the  result  of  applying  the  same  double  optimization 
procedure  to  the  Poisson  noise  image  sets,  with  the  result  that  the  .tss  error 
between  thr  predicted  and  observed  quality  is  minimized  at  3.0  Inches  rms  for 
an  a of  1125;  it  should  be  noted  that  this  error  is  somewhat  less  than  the 
3.5-inch  rms  value  found  using  the  eye  MTF  independent  model. 

4.3.5  Conclusions  Concerning  Fit  to  Complete  Summary  Measure  Model 

As  shovn  above,  the  use  of  the  complete  summary  measure  model  yields 
some  improvement  over  the  performance  of  the  eye  MTF  independent  model  based 
on  a root  mean  square  error  criterion.  However,  the  use  of  the  complete 
model  has  resulted  in  a significant  reduction  in  the  average  fractional  er- 
ror, as  demonstrated  by  th?  reduction  in  the  scatter  of  the  results  for  small 
values  of  quality. 

If  the  summary  measure  model  is  expressed  as 

Q - DQRF  * Qq  (37) 

where 

Q - predicted  quality 
Qq  - quality  based  on  OTF  only  (no  noise) 

DQRF  ■ dimensionless  quality  reduction  factor 

it  can  be  seen  that  the  dimensionless  quality  reduction  factor  (DQRF)  has  a 
much  smaller  error  for  the  case  of  low  quality  estimates  than  dees  the  cor- 
responding DQRF  for  the  simple  model. 
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4.4  USE  OF  MODEL  FOR  SYSTEM  PERFORMANCE  ANALYSIS 


; 


In  the  previous  paragraphs,  it  was  shown  that  the  simple  eye  MTF  indepen- 
dent summary  measure  model  was  sufficient  to  predict  image  quality  to  within 
an  rms  error  of  less  than  4 inches;  in  this  paragraph,  this  model  will  be  used 
in  a sample  calculation  to  investigate  the  tradeoff  between  image  quality,  rms 
wavefront  error,  and  exposure  time. 


By  use  of  Equation  (29) , the  resultant  image  quality  for  a given  set  of 
parameters  can  be  derived;  conversely  the  exposure  time  for  a specified  image 
quality  can  be  derived  by  rewriting  Equation  (29)  as  follows: 


t * 

e 


4 2 

D S 
o 


2 M2 

a a Nt 
r I 


(38) 


where 


N ■ noise  standard  deviation  increase  factor  resulting  from 
post-compensation  restoration  filter 

S * normalized  Strehl  ratio  (normalized  to  1.0  for  diffraction- 
limited  system) 

Figure  29  illustrates  the  use  of  Equation  (38)  to  derive  a set  of  three 
isoquality  curves,  for  an  arbitrary  object,  for  long  time  exposure  imaging. 

In  Figure  29,  each  point  on  an  isoquality  contour  represents  a combination 
of  exposure  time  and  pre-compensation  system  residual  error  that,  when  com- 
bined with  a post-processing  operation  that  yielded  a Strehl  ratio  of  0.8, 
results  in  a system  image  quality  of  the  denoted  value.  It  should  be  noted 
that  a post-compensation  Strehl  ratio  of  0.8  corresponds  to  a quality  Q of  54 
inches;  i.e.,  noiseless  imagery  would  have  the  subjective  quality  of  images 
obtained  through  a diffraction-limited  54-inch  aperture.  Also,  it  should  be 
noted  that  as  the  residual  wavefront  error  increases,  the  exposure  time  re- 
quired to  yield  the  same  resultant  image  quality  also  increases. 
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CONCLUSIONS,  RECOMMENDATIONS  AND  SUMMARY 


5.1  CONCLUSIONS  CONCERNING  FIT  TO  SIMPLE  MODEL 

As  a result  of  the  application  of  the  simple  summary  measure  model  to 
the  total  number  of  106  images  generated  on  the  study,  the  simple  model  was 
found  to  predict  quite  well  the  observed  quality  of  the  images,  for  a wide 
range  of  noise  and  OTF  conditions. 

The  scale  parameter  a was  found  to  be  1800  for  the  66  ’"oisson  noise 
images,  and  resulted  in  a 3.5-inch  rms  difference  between  the  subjective 
quality  (as  determined  by  the  viewing  and  interpretation  experiments) , and 
the  quality  predicted  by  the  model.  For  the  30  Gaussian  images,  a value  of 
a - 3400  was  found  to  yield  a fit  with  3.3-inch  rms  error;  for  the  10  noise- 
less images,  a fit  of  2.9  inches  rms  was  found. 

The  simple  model  was  found  to  have  a larger  percentage  error  in  the 
region  of  poor  image  quality,  and  was  found  to  perform  far  better  on  a 
percentage  basis  for  the  case  of  good  image  quality,  which  is  precisely  the 
region  where  the  tool  will  be  most  needed  -i.e.,  for  tradeoffs  concerning 
the  performance  of  a compensated  imaging  system  that  is  producing  images 
not  too  far  from  the  quality  of  diffraction  limited  images. 

Thus,  it  is  concluded  that  the  simple  model  will  serve  well  as  a tool 
for  estimating  the  quality  of  satellite  imagery,  for  the  type  of  generic 
degradations  simulated  in  this  study. 


5.2  CONCLUSIONS  CONCERNING  FIT  TO  COMPLETE  MODEL 

The  application  of  the  more  complete  summary  measure  model,  including 
the  eye  magnification  optimization  step,  resulted  in  a moderate  reduction 
in  the  rms  error  of  the  fit  of  the  model  to  the  results  of  the  V and  I ex- 
eriments.  For  the  Poisson  noise,  a value  of  a of  1125  yielded  an  rms  error 


of  3.0  inches,  while  for  the  Gaussian  images,  a value  of  a of  2000  yielded 
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an  rms  error  of  3.5  inches.  The  fit  of  the  noiseless  images  was  reduced 
from  2.9  inches  rms,  to  2.0  inches  rms. 

However,  the  complete  model  performed  significantly  better  on  a per- 
centage basis,  as  indicated  by  the  reduction  of  the  scatter  of  the  data 
shown  in  Figures  27,  28,  and  29,  in  the  region  of  poor  image  quality. 

Thus  it  is  concluded  that  the  complete  model,  including  the  eye  transfer 
function,  does  provide  a more  complete  description  of  the  relationship 
between  the  perceived  image  quality  (as  determined  by  human  observers) , and 
various  descriptions  of  system  transfer  function  and  noise  performance. 


5 . 3 RECOMMENDATIONS 

As  a result  of  the  study  described  in  this  report,  two  recommendations 
are  made,  as  follows: 

(1)  Repeat  some  of  the  experiments  performed  during  the  study  for 
several  different  satellite  Images. 

(2)  Examine  the  question  of  optimum  post-processing  by  generating 
images  that  have  been  post-compensated  by  Weiner  filters, 
geometric  mean  filters,  and  any  other  algorithms  that  appear 
to  have  applications  to  the  post -compensation  function. 

5.4  SUMMARY 

In  summary,  it  is  conceded  that  a useful  model  exists  relating  the 
subjective  image  quality  of  satellite  imagery  to  various  quantitative 
measures  of  system  performance;  the  quantitative  measures  that  are  required 
include  the  system  optical  transfer  function,  and  the  spectrum  and  strength 
of  the  noise  degrading  the  Imagery. 

With  this  model  in  hand,  it  is  now  possible  to  perform  system  tradeoffs 
to  determine  the  optimum  split  between  the  performance  of  the  pre-comp- 
ensation function  and  the  post-compensation  function  for  a compensated 
imaging  system;  in  addition,  it  is  possible  to  perform  tradeoffs  concerning 
the  form  and  nature  of  the  post-prpeessing  algorithms  to  be  used  on  the  post 

.compensation  function. 
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SECTION  A- I 


INTRODUCTION 

This  appendix  describes  a technique  for  generating  random  wavefronts 
that  simulate  the  statistics  of  an  atmospherically  distorted  wavefront.  The 
atmosphere  is  described  by  a Kolmogoroff  spectrum  and  the  wavefront  is  treated 
as  an  expansion  in  Zernike  polynomials. 

Section  A-II  describes  the  mathematical  preliminaries.  Section  A-III 
states  the  procedure  for  generating  wavefronts,  and  Section  A-IV  contains  the 
evaluation  of  an  integral  used  in  Section  A-II. 
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SECTION  A-II 


ATMOSPHERIC  STATISTICS  FOR  A MODAL  REPRESENTATION 


The  aberration  function  $ will  be  expanded  in  terms  of  a complete  set  of 


functions  { f ^ ( pQ ) } . 


4>(p0)  - Z a f (p6) 
j J J 


(A-l) 


The  functions  f are  orthonormal  so  that 

1 2tt 


| J pdp  d0  f*  (P0)  fj,  (P0)  - fijj, 


(A- 2a) 


Jd2r  f*  (r)  f j , (r)  - 


(A- 2b) 


The  mean  square  error  E produced  by  such  an  expansion  Is  defined  as 


- J d2r  |$(r)  - $n(r)J 


(A-3) 


where  4>  (p0)  is  the  nfc  order  polynomial  expansion  in  Equation  (A-l).  Sub- 
n 

stituting  Equation  (A-l)  info  Equation  (A-3)  for  #n  and  finding  the  set  of 
(dj)  that  minimizes  E yields 


E - f d2r  4>2(r)  - Z |a  |‘ 

j-1  J 


(A-4) 


a . - J d2r  4>(r)  fj(r) 


(A-5) 
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The  aberration  function  is  taken  to  be  that  associated  with  an  atmospherically 
distorted  random  wavefront  so  that  an  average  mean  square  error  can  be  defined 
as 

<E  > - J*  d2r  <$2(r)  > - E < |c.j|2>-A  <A_6) 

Over  the  aperture  the  random  wavefront  is  assumed  stationary  so  that  Equation 
(A-6)  reduces  to 

A - to2  - E < |a^|2>  , (A-7) 


where  a2  is  the  variance  of  the  wave  fluctuations  and  can  be  defined  in 
terms  of  the  power  spectral  density  W(k)  of  the  fluctuations  as 


2 

0* 


! 


d2k 

(27T)2 


W(k) 


(A-8) 


where  k = (k  ,k  ) 
x y 

The  nean  square  value  of  the  coefficients  {c^}  can  now  be  expressed  in 
terms  of  W(k) . The  Fourier  transform  of  the  functions  f^(r)  is  defined  as 


F^(k)  = J*  d2r  f j <r)  e 


-ik'r 


(A-9) 


where  r = (x,y) 


with  inverse  transform 


f (D  - F (k)eik-r 

J (2tt r 3 


(A-10) 


From  Equation  (A-f)  , the  mean  square  value  of  the  coefficients  is 


<|aj|2>  = JJ*  d2r'd2r  <$(r)  $(r')  > f j (r)  f^r’) 


(\-ll) 
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Substituting  Equation  (A-10)  into  Equation  (A-ll)  and  making  the  transforma- 
tion, 


r'  - r + t 


yields 

2 

< ja  |2  > - | F , (k)  | 2 J*d2x  <4*(r)  *(r-K)  > ei(k’T)  (A-12) 

3 (2  it)  ^ 3 

Equation  (A-12)  is  thus  easily  expressed  in  terms  of  the  fluctuation  spectrum 
as 

<|aj2>  - f-^r  |F  (k)  | 2 W(k)  (A-13) 

J J (2ir)Z  2 


Generalizing  Equation  (A-13) , the  quantity  < > can  be  written  as 


,> 


.2. 

(*-—  F*  (k)  W(k)  F ,(k) 
J (2n)2  3 3 


(A-14) 


Equation  (A-14)  defines  a matrix  represtntation  of  the  spectrum  W(k)  which  is 
denoted  by 


(A-15) 


One  further  item  should  be  mentioned  In  connection  with  a generalized 
K>dal  expansion,  and  that  is  closure.  Equation  (A-7)  describes  the  error  re 
suiting  from  taking  a finite  number  of  terms  in  the  expansion.  The  quantity 
A must  always  be  positive  and  as  such  provides  a good  measure  for  deciding 
how  efficiently  the  random  wavefront  hes  been  generated.  An  alternative 
form  for  Equation  (A-7)  is 

A - W(k)  - Z |F.  (k)  | 2 W(k) 

(2it)2  j (2tt)Z  3 


(A-16) 
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Zernike  Polynomial  Expansion 

In  the  following  the  orthonormal  set  of  functions  fj(p0)  are  taken  to 
be  Zernike  Polynomials,  Z^(pG)  where 


tm 


Z,(p0)  = f (pa) 

J n 


Vn  + 1 „±m  , N ±im0 

— Rn  (p)  6 


(0  < pi  ) 


(A-17) 


The  first  few  polynomials  are  shown  in  Table  A-l  with  Z^  = fj/VW.  The 
Fourier  transform  of  Equation  (A-17)  is 


j (k) 

Q (k)  - F“  (k^g)  = 2im^Tl(n+l)  (-l)(n-m)/2  eimk9 


(A-18) 


The  spectrum  of  turbulence  is  taken  to  be  a Kolmogoroff  spectrum,  which  can 
be  written  as 


„ \5/3  11/1 
W(k)  = 19.29  ( j-  1 k 11/3 

o 


(A-l 9) 


Equation  (A-19)  is  different  from  the  usual  Kolmogoroff  spectrum  in  that  the 
radius,  a,  of  the  telescope  aperture  is  included.  This  definition  is  con- 
venient because  it  allows  all  equations  and  computer  calculations  to  be  per- 
formed for  a unit  radius  aperture.  Thus,  the  spatial  frequency  k has  units 
of  radians/unit  radius. 

The  matrix  representation  defined  in  Equation  (A-15)  yields  a Zernike 
representation  of  the  Kolmogoroff  spectrum  written  as 


ajj'  = 2 


19.29 


(%)*'*  (n+1)(n’+1) 


1/2 


(-1) 


(n+n'-2m)/2 


; 


dk 


Wk)  Jn’+l(k) 
, 14/3 


(A-20) 


Because  of  the  relationship  between  the  j indices  and  the  n,m  indices  implied 
in  Table  A-l,  the  matrix  elements  ajj»  with  even-odd  subscripts  are  zero 
(i.e.,  s * 0,  etc.).  This  is  a consequence  of  the  circular  symmetry 
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in  W(k) . The  inodes  associated  with  odd  j have  their  angular  dependence  de- 
scribed by  sin  m 0 while  those  of  even  j are  described  by  cos  m 0.  Clearly, 
the  angular  Integrations  with  mixed  sines  and  cosines  are  zero.  The  integral 
in  Equation  (A-20)  is  a tabulated  integral  and  yields  a number  for  all  the 
matrix  elements  a^,  except  a^  (n*n*  ■ 0> , the  piston  coefficient.  Since 
the  total  variance  of  the  Kolmogoroff  spectrum  is  infinite,  it  is  not  sur- 
prising to  find  an  infinite  piston  coefficient.  Leaving  the  integral  evalua- 
tions to  Section  IV,  it  is  found  that  the  closure  relation.  Equation  (A-16) 
for  the  piston  corrected  spectrum  can  be  written  for  the  Zernike  representa- 
tion of  the  Kolmogoroff  spectrum  as 


dim 

k -K) 
o 


it  J H(k)  - J |Q1(k)  [ 2 W(k) 


kdk 

2tt 


(A-21) 
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or  that 


— = 1.03 
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it) 
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(rad2) 


(A-22) 


which  is  the  mean-square  piston  corrected  error  first  derived  by  Fried  . 
Thus,  infinite  variance  of  the  spectrum  is  associated  only  with  the  piston 
and  is  of  no  serious  physical  consequence. 


With  the  value  of  the  integral  in  Equation  (A-20)  given  in  Section  IV 
the  exact  mean-square  error  resulting  from  using  N modes  can  be  evaluated. 
In  the  limit  of  large  N,  an  asymptotic  formula  has  been  extracted, 


— - 0.2944  N~0,866  (D/r  )5/3  (rad2) 

TT  O' 

*D.L.  Fried,  J.  Opt.  Soc.  Amer.,  55  1427  (1965). 


(A— 23) 
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The  matrix  elements  given  by  Equation  (A-20)  are  uaed  in  the  next  section  to 
generate  random  wavefronts. 
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SECTION  A-III 


GENERATING  RANDOM  WAVEFRONTS 


In  this  section  the  method  of  generating  random  wavefronts  is  discussed. 
The  heart  of  the  method  rests  upon  the  matrix  A described  in  Section  A-II. 
Because  the  Kolmogoroff  spectrum  is  a real  symmetric  spectrum  the  matrix  A 
is  a real  symmetric  matrix  and  as  such  can  be  factored  into 


A * L L 


(A-24) 


where  L is  a lower  triangular  matrix 


L is  the  transpose  of  L 


The  Cholesky2  decomposition  method  of  factoring  A was  used  in  all  random  wave- 


front  simulations. 


The  object  of  a wavefront  simulation  is  to  generate  random  numbers  {a^} 


such  that 


<a,  > =0,  <aV,  > 


^ ' jj’ 


(A-25) 


This  can  be  achieved  if 


E q 


(A-26) 


<qt>  =o,  <qiqj> 


(A-27) 


The  desired  random  wavefront  is  then  from  Equation  (A-l) 


4>(p0)  = £ o,  f (p0) 

J J 


(A-28) 


2See,  e.g.,  J.H.  Wilkinson,  "The  Algebraic  Eigenvalue  Problem",  Oxford, 
Fairlawn,  New  Jersey  (1965). 
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A computer  generated  random  wavefront  then  consists  of  the  procedures,  as 
shown  in  Figure  A-l. 

(1)  Calculate  the  matrix  A. 

mtm 

(2)  Factor  the  matrix  A into  L L. 

(3)  Generate  the  random  numbers  q^. 

(4)  Calculate  random  coefficients  from  Equation  (A-26) . 

(5)  Calculate  random  wavefront  from  Equation  (A-28) . 
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MONTE  CARLO  SIMULATIONS  OF 


Random  Wavefront  Simulation 


r 


SECTION  A-IV 


EVALUATION  OF  THE  INTEGRAL  I , 

nn 


In  this  section  the  Integral  in  Equation  (A-20)  I , is  evaluated. 


‘nn-  ■ J 


dk 


Wk>  Jn’+l(k) 

0473 


(A-29) 


This  Integral  is  tabulated  in  Most  standard  integral  table  handbooks  of  Bessel 
function  Integrals.  ^ 


r««/3>  -2u/?— ) 


/ -n+n'  + 14/3  + 1 \ p / 

f n-n'  + 14/3  + 1\  r i 

f'n+n'  +14/3  + 3 \ 
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l 2 ) r 1 

l 2 / 

for  n,n'  i 0 

where  T(x)  is  the  gaasa  function. 

The  piston  Integral  I f as  Indicated  in  the  text,  needs  special  attention. 
Let  us  write 


oo 
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| Jx(k) |2  dk 
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(A-31) 


Equation  (A-31)  can  be  Integrated  by  using  the  Identity 
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Equation  (A-31)  then  can  be  written  as 


roo  * if  J ' dr  ’<*> 


2tt  00 

J df)  j 

o k 


, -8/3  -ikr  cos0 
dk  k e 


The  k integration  is  recognized  as  an  incomplete  gamma  function  so  that 


.-8/3  -ikr  cos0 
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The  first  term  on  the  right  of  Equation  (A-35)  can  be  written  as 


(ir  cos0)5/3  r(-5/3)  i3(r  ~os0)3  J t4/3_1  e"lrt  COS0 


Substituting  Equations  (A-35)  and  (A-36)  into  Equation  (A-34)  and  performing 
the  angular  integration  yields 

00 

I .-if  rdtg(r)  M k -5/3.271.  1/3  5/3  fs  [(7/6)1 
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The  Integrals  over  r are  easily  expressed  in  terns  of  the  beta  function 
B(x,y)  by  means  of 

1 »“  |^co.-lv-v  £77  dv - i [-Ij. - -ij j,(^.  i)  ( 


(A- 38) 


The  result  Is 


I . _L  k -5/3  3888  T(7/6)  1(7/3.  1/2] 
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other  terms  (A- 39) 


The  other  terms  in  Equation  (A-39)  vanish  when  Kq  Is  aero,  and  the  constant 
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Figure  B-2.  Poisson  Noise  with  Various  Diffraction- 
Limited  Apertures  (Set  B) 
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Figure  B-3.  Gaussian  Noise  with  Various  Diffraction- 
Limited  Apertures  (Set  H) 
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Figure  B- 




. Short  Time  Exposure  OTF  and 
Poisson  Noise  (Set  C) 


90 


129  130  131 


w 


150 


151 


15'. 


Figure  B- 


. Short  Time  Exposure  OTF  and  Gaussian  Noise 
and  Post-Processing  (Exact)  (Set  K) 
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Figure  B-10.  Short  Time  Exposure  OTF  and  Poisson  Noise  and 
Post-Proc easing  (Approximate)  (Set  J) 
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APPENDIX  C 


SET  BY  SET  SUMMARY  OF  RESULTS 
OF  FITTING  SIMPLE  MODEL  TO  SIMULATION 
STUDY  IMAGES 
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The  eye  modulation  transfer  function  used  in  the  full  summary  measure 

1 2 

model  is  *n  analytical  approximation  of  published  data  ’ given  by 


T 

eye 


00 


U0.8k 
1 + 0.025  k4 


where  k = spatial  frequency  (cy/mm) 
Figure  D-l  illustrates  this  function. 


Figure  D-l.  Eye  Transfer  Function 


1-0.  Bryngdahl,  Optika  Acta,  _12,  Jan.  65 
2 - J.J.  DePalma,  E.M.  Lowry,  J0SA,  52,  Mar.  62,  p 328 
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APPENDIX  E 

ESTIMATIVE  FILTER  RESULTS 
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An  important  consideration  is  the  application  of  the 
Summary  Image  Quality  (SIQ)  model  to  generalized  restoration 
filters.  In  order  to  qualify  the  performance  of  the  SIQ  model 
two  generalized  filters  were  used  to  restore  long  time  exposure 
images.  These  two  filters  were  the  Wiener  and  Geometric  Mean 
filter. 


The  Wiener  filter  is  the  best  generalized  filter  in  a 
mean-square-error  sense.1  This  filter  takes  the  form: 


F(k) 


T*fi(k) 


I Tj.e  (k>  I 2 + Vk)/Vk) 


E-l 


where 

$n  (k ) is  the  estimate  of  the  noise  power  spectrum 
*n(k)  is  the  estimate  of  the  object  power  spectrum. 


The  Geometric  Mean  filter  is  a geometrical  mean  between 
the  inverse  filter  and  the  Wiener  filter.  This  filter  as  sug- 
gested by  Stockham,  Cole,2  and  Cannor3  is,  in  its  general  form: 


F (k) 


T*  (k) 

le 


♦ *«<k>/*  <k) 

it  0 


l-s 


E-2 


*C.  w.  Helstrom,  "Image  Restoration  by  the  Method  of  Least 
Squares,"  JOSA,  Vol.  57,  No.  3,  pp  297-303,  (Mar  1967). 

2E.  R.  Cole,  "The  Removal  of  Unknown  Image  Blurs  by  Homomorphic 
Filtering,"  Ph.  D.  Dissertation,  Dept,  of  E.  E.,  Univ.  of  Utah, 
Salt  Lake  City,  (Jun  1973). 

*T.  M.  Cannor,  "Digital  Image  Deblurring  by  Nonlinear  Homomor- 
phic Filtering,"  Ph.  D.  Dissertation,  Dept,  of  E.  E.,  Univ.  of 
Utah,  Salt  Lake  City,  (Aug  1974). 
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where  0<S^1.  For  the  isotropic  phaseless  filter,  II-H*  &nd  S=l/2, 
equation  E-2  reduces  to: 


F (k)  = 


* 00 

o 


♦ (k)  |Tle(k)  !2  + *n(k) 


1/2 


E-3 


again  for  the  lonq  time  exposure  case.  Figures  E-l  and  E-2  re- 
present various  degrees  of  wavefront  error  and  noise  variance 
for  the  Wiener  and  Geometric  Mean  cases  respectively. 

When  the  data  from  V and  I experiments  were  reduced, 
agreement  of  the  SIQ  model  to  observed  image  quality  was  within 
2.5  inches  standard  deviation. 


114 


Figure  E-l  SET  L.  Long  Time  Exposure  Atmospheric  OTF  and  Poisson 
Noise  and  Vfiener  Filter  Restoration  to  Strehl 
Ratio  ■ 0.8. 
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Figure  E-2  SET  M.  Long  Time  Exposure  Atmospheric  OTF  and  Poisson 
Noise  and  Geometric  Mean  Filter  Restoration  Strehl 
Ratio  * 0.8. 
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GLOSSARY 


a - photoelectrons  per  pixel  per  unit  reflectivity 

A ■ coefficient  of  nth  Zernike  polynomial 
n 


c 

C(r) 


dimensionless  scale  factor 

two-dimensional  autocorrelation  function  of  wavefront 


DQRF 

D 

mse 

D 

o 

D(r) 


dimensionless  quality  reduction  factor 
mean  square  error  image  quality  criterion 
aperture  diameter  (inches  or  meters) 
phase  structure  function 


G ■ edge  response 

E — mmriwmm  value  of  filter  function 
max 


F 

F (k) 
8 

F 

Fe(k) 

f 


focal  length  (meters) 

approximate  transfer  function  restoration  filter 
fidelity  criterion 

exact  transfer  function  restoration  filter 
relative  fidelity  criterion 


i(r) 


image 


J^(2uRk) 


k 

k 


nyq 

*(x) 


Bessel  function 

(k  ,k  ),  two-dimensional  spatial  frequency  (cycles/milliaeter) 
x y 

normalized  spatial  frequency 
Nyquist  frequency 

line  spread  function 


m 


”0 


viewing  magnification 
optimum  viewing  magnification 
pre-noise  image  modulation  factor 
post-noise  image  modulation  factor 
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GLOSSARY  (Continued) 


n(x) 

NI 

N(k) 

NRE 


noise  random  variable 

additive  Gaussian  noise  standard  deviation  increase  factor 
noise  power  spectral  density 
normalized  relative  energy 


o(r) 


object 


P 

Q 

^ioc 

q„ 


point  spread  function 

image  quality  (inches) 
image-object  correlation  criterion 
image  quality  for  no  noise 


s 

R 


(u,v)  two-dimensional  coordinate  on  aperture  plane  (meters) 
coherence  length 
reflectivity  of  satellite 
aperture  radius  (inches) 


Strehl  ratio 


T0(k) 

Tdt<k> 

Teye(k) 

Tls 

Tka(k> 

Vk> 

T <k) 
o 

T (k) 
se 


Ug(x) 


Vk) 


relative  image  sharpness 

exposure  time 

display  transfer  function 

desired  diffraction-limited  transfer  function 

eye  modulation  transfer  function 

image  sharpness  criterion 

known  degrading  transfer  function 

long  time  exposure  atmospheric  transfer  function 

diffraction-limited  transfer  function 

short  time  exposure  atmospheric  transfer  function 

value  of  sampled  image  point  (pixel)  (photoelectrons) 

2 9 

power  spectrum  of  atmospheric  turbulence  (rad  /(cy/m)  ) 
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GLOSSARY  (Continued) 


Z (u,v)  ■ nth  Zernike  polynomial 

n 

a ■ noise  scaling  factor 
6 - filter  parameter 

e ■ mean  square  wavefront  error  (radian  ) 

H - fidelity  criterion  parameter 

0 - pixel  angular  subtense  (radians) 

A - wavelength  (micrometers) 

o - image  noise  (photoelectrons) 

O * image  noise  (reflectivity) 

■ average  image  noise 

$(r)  ■ wavefront  phase  (radians) 

*(k)  ■ object  power  spectrum 

ft  - Nyquist  angular  sampling  rate  (samples/radian) 
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MISSION 

of 

Rome  Air  Development  Center 


r adc  is  the  principal  AFSC  organization  cnargea 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences , intelli- 
gence, coaaoand,  control  and  communications  technology, 
products  and  services  oriented  to  the  needs  of  the  USAF. 
Primary  RADC  mission  areas  are  communications , electro- 
magnetic guidance  and  control,  surveillance  of  ground 
and  aerospace  objects,  Intelligence  data  collection  and 
handling,  information  system  technology,  and  electronic 
reliability > maintainability  and  compatibility.  RADC 
has  missio, : responsibility  as  assigned  by  AFSC  for  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence,  mapping,  charting,  command 
control  and  coamtunl ca tions  areas. 


